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Abstract. The development of germanium Compton 
telescopes for nuclear 7-ray astrophysics (^0.2-20 MeV) 
requires new event reconstruction techniques to accu- 
rately determine the initial direction and energy of photon 
events, as well as to consistently reject background events. 
This paper describes techniques for event reconstruction, 
accounting for realistic instrument/detector performance 
and uncertainties. An especially important technique is 
Compton Kinematic Discrimination, which allows proper 
interaction ordering and background rejection with high 
probabilities. The use of these techniques are crucial for 
the realistic evaluation of the performance and sensitivity 
of any germanium Compton telescope configuration. 
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1. Compton telescopes for 7-ray astrophysics 

Looking beyond the INTErnational Gamma-Ray Astro- 
physics Laboratory (INTEGRAL), the next generation 
soft 7-ray (~0.2-20 MeV) observatory will require high 
angular and spectral resolution imaging to significantly 
improve sensitivity to astrophysical sources of nuclear 
line emission. Building upon the su ccess of COMP- 
TEL/CGRO flSchonfelder et al. 1993|) , and the high 
s pectral resolution o f the upcoming SP I/INTEGRAL 
( |Vcdrcnne ct al. 1998], Lichti ct al. 1996 ) , a number 



of research ers ( Johnson et al. 199q , Jean et al. 1996 , 
Boggs 1998 ) have discussed the merits of a high spec- 
tral/angular resolution germanium Compton telescope 
(GCT); the ability to achieve high sensitivity to point 
sources while maintaining a large field-of-view make a 
high resolution Compton telescope an attractive option 
for the next soft 7-ray observatory. 

The development of Compton telescopes began in 
the 1970's, with work done at the Max Planck Institut 



side (Herzo et al. 1975), and the University of New Hamp- 
shire ^Lockwoo^^^lrT97^) , culminating in the design 
and flight of COMPTEL/CGRO. These historical Comp- 
ton telescopes consist of two scintillation detector planes - 
a low atomic number 'converter' and a high atomic num- 
ber 'absorber.' The model interaction of a Compton tele- 
scope is a single Compton scatter in the converter plane, 
followed by photoelectric absorption of the scattered pho- 
ton in the absorber. By measuring the position and energy 
of the interactions, the event can be reconstructed to de- 
termine the initial photon direction to within an annulus 
on the sky. 

A handful of groups are actively developing imaging 
germanium detectors (GeDs) partly in anticipation of a 
GCT ( |Luke et al. 1994] |Kroeger et al. 1996| ). The goal of 
these researchers is to develop large area detectors with 
(sub) millimeter spatial resolution, while maintaining the 
high spectral resolution (E/SE ~ 500 at 1 MeV) charac- 
teristic of GeDs. The use of high spectral/spatial resolu- 
tion GeDs as converter and absorber planes would signif- 
icantly improve the performance of a Compton telescope, 
but will add a number of complications to the event recon- 
struction. Most significantly, with the moderate atomic 
number (Z = 32) of germanium, photons will predomi- 
nantly undergo multiple Compton scatters before being 
photoabsorbed in the instrument. Furthermore, with in- 
teraction timing capabilities of ~10 ns, the interaction 
order will not be determined unambiguously by timing 
alone. Compton Kinematic Discrimination (CKD) is pro- 
posed here to overcome these complications, an extension 
of a method first discussed in context of liquid xenon 



( Bchonfelder et al. 1973 ), University of California, River- 
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time projection chambers (Aprilc ct al. 1993). The abil- 
ity of this technique to allow proper event reconstruction 
is investigated in detail. 

Due to their relatively low efficiency (typically ~ 1%), 
Compton telescopes rely on efficient background suppres- 
sion to maintain their sensitivity. In addition to inter- 
action ordering, techniques are presented using CKD, in 
combination with other tests and restrictions, to suppress 
the dominant background components. 
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Fig. 1. Example Compton telescope. If a photon under- 
goes one or more Compton scatters in the instrument and 
then is photoelectrically absorbed, then by using the po- 
sitions (n, . . . , t\at) and energy deposits (E±, . . . , En) of 
the interactions, the initial direction of the photon can 
be determined by the Compton scatter formula to within 
an annulus on the sky, <j>i. The width of this annulus is 
determined by the uncertainties in both the interaction 
locations and energy deposits. 



The goal of this work is to outline a complete set of 
event reconstruction techniques for GCTs, taking into ac- 
count realistic detector/instrument performance and un- 
certainties. Examples of the techniques are presented for a 
GCT configuration outlined in Appendix A; however, full 
analysis of this configuration will be presented in a second 
paper dedicated to the optimization and performance of 
several GCT configurations. The full analysis of a GCT 
configuration is complicated, requiring a detailed study of 
the tradeoffs between efficiency, angular and spectral res- 
olution; therefore, this paper focuses only on the detailed 
discussion of the event reconstruction techniques which 
will be used in future work dedicated to analyzing GCT 
performance. 



2. Principles of Compton imaging 

The principle of Compton imaging of 7-ray photons is 
illustrated in Figure 1. (See von Ballmoos, Diehl and 
Schonfelder (1989) for an excellent review of historical 
Compton telescope configurations.) An incoming photon 
of energy E and direction p undergoes a Compton scatter 
at an angle <j>\ at the position n within a detector, creat- 
ing a recoil electron of energy of E\ which is quickly ab- 
sorbed and measured by the detector itself. The scattered 
photon then deposits the rest of its energy in the instru- 
ment in a series of one or more interactions of energies Ei 
at the positions r i: until eventually photoabsorbed. Here 
the total photon energy after each scatter i, normalized to 
the electron mass, is defined as 



Wi = 



1 



m e c z 



N 
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j=i+l 



Ei 
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where W a — E/m e c 2 , and N is the total number of inter- 
actions. The initial photon direction is related to scatter 
direction vector = r-i — T\ (r\ after normalization), and 
the scattered photon energies Wj by the Compton formula 

: = ' - TTt 777" ■ (2) 
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Given the measured scatter direction and the angle 
cos</»i implied from the energy depositions, the equation 
for p is not unique (if the electron recoil direction could 
be measured, it would solve this ambiguity); therefore, 
the initial direction of the photon cannot be determined 
directly, but it can be limited to an annulus of directions 
p' which satisfy the equation 



r 1 ■ p — cos tpi 



(3) 



There are two uncertainties in determining the event 
annulus: the uncertainty in <f>i due to the finite energy 
resolution of the detectors, here labelled and the 

uncertainty in r[ determined by the spatial resolution of 
the detectors. Both of these uncertainties add to deter- 
mine the uncertainty (effective width) of the event annulus 
S(f>i. From Equation 2, the derivation of 8<j)\ t E is straight- 
forward and yields 



where, 
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In order to simplify the analysis, it is convenient to 
transform the uncertainty Sr^ into an effective uncertainty 
in (pi, defined as 5</>i ir , such that 
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Fig. 2. Illustration of how uncertainties in the scattered 
photon direction r' can be translated into effective uncer- 
tainties in the scattering angle <fi. 



The angular resolution 8<f>i tT is the effective 'wiggle' of 
r'i around its measured direction due to the uncertain- 
ties in the spatial measurements. The spatial uncertainties 

are defined as 5x[ — yj Sxf + Sxf +1 , 8y[ = y 1 Syf + Syf +1 , 



Szl = yj Szf + Szf +1 . It is simplest to analyze the situation 
for each axis separately as shown in Figure 2. The uncer- 
tainty in the direction of due to the uncertainty 5x[ is 
given by 



tan(^) = (M)^l-(r^ 



(7) 



Likewise for the other axis, 



which combine to yield the total uncertainty S(j)i r given 

by 



+ . 



(8) 



For detectors with a given energy resolution, in order 
to optimize the performance of a Compton telescope one 
would require that < 54>i,E in the energy range of 

interest. To first order, this implies that the spatial reso- 
lution in relation to the scale size of the instrument must 
be comparable to or less than the energy resolution, i.e. 
Sri/r', < SE[/E[. 



3. Complications of germanium Compton 
telescopes 

Finite detector thresholds, energy resolutions, and spatial 
resolutions produce systematic biases in the imaging ca- 
pabilities of Compton telescopes. These limitations have 
been discussed in detail elsewhere in context of two-layer, 
low-Z converter and high-Z absorber, scintillation detec- 
tor designs (von Ballmoos ct al. 1989), and the conclu- 
sions can be directly applied to GCTs. However, GCT de- 
signs will introduce additional complications which signifi- 
cantly alter the event reconstruction techniques. Historical 
Compton telescope configurations make two assumptions 
about the events which do not generally hold in GCTs: (i) 
the events are a single Compton scatter in the converter, 
followed by photoelectric absorption in the absorber, and 
(ii) the time-of-flight (TOF) between the photon interac- 
tions is measured to determine their order. 

The distributions of number and type of interaction 
sites in a GCT for normally incident, fully-absorbed pho- 
tons ranging from 0.2-10 MeV are shown in Figure 3, for 
the instrument configuration discussed in Appendix A. 
Here we distinguish three event types: a single photoelec- 
tric absorption, one or more Compton scatters followed 
by a single photoabsorption, and one or more pair pro- 
ductions. Compton scatters followed by pair production 
could potentially be reconstructed; however, here we in- 
clude these events with other pair productions. These dis- 
tributions account for the finite spatial resolution of the 
detectors, so that interactions occurring too closely to- 
gether are not resolved. From these distributions it is clear 
that events with ~8 or more interaction sites can be imme- 
diately rejected as probable pair production events, with 
little effect on the Compton photopeak efficiency. For in- 
cident photon energies above 0.5 MeV, 3-7 interaction site 
Compton scatter events dominate the photopeak. 

To accurately reconstruct a Compton scatter event, 
the first and second interaction sites must be spatially 
resolved, and their order correctly determined. The need 
to determine the proper ordering of three or more (3+) 
interaction sites is complicated by the timing capabil- 
ities of GeDs. In the scintillation detectors of COMP- 
TEL/CGRO the interaction timing can be performed to 
~0.25 ns (Schonfeldcr et al. 1993), which is adequate to 
determine the TOF between two interactions in the sep- 
arate detector planes. With the slower rise time of GeDs 
one can reasonably expect event timing to ~10 ns, which 
is inadequate for TOF measurement in reasonably-sized 
instruments. While Pulse Shape Discrimination methods 
have been proposed to push the interaction timing in GeDs 
to ~1 ns ( [Boggs 199E ), even this timing would be unreli- 
able for determining TOF among three or more interac- 
tion sites. A method of reliably determining the photon 
interaction order without timing information must be de- 
veloped. 
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Fig. 3. Statistical distributions of the number of interac- 
tion sites for fully-absorbed photons for the instrument 
configuration presented in Appendix A. These distribu- 
tions take into account the finite spatial distribution of 
the detectors, combining interactions that cannot be spa- 
tially resolved. Events are divided into photoelectric ab- 
sorptions (solid white), pair productions (solid black), and 
the desired Compton scatter (s) followed by a single pho- 
toabsorption (striped). Events with 12 or more sites were 
combined in a single bin. 



4. Multiple Compton scatter events: Compton 
kinematic discrimination 

One method has been suggested to overcome these com- 
plications in the context of liquid xenon time projection 
chambers (Aprile et al. 1993). Here, this method is for- 
malized as Compton Kinematic Discrimination (CKD) 
and examined in more detail. This technique allows the 
order of the photon interactions to be determined with 
high probability, as well as providing the basis of a pow- 
erful tool for background suppression in GCTs. 

CKD takes advantage of redundant measurement in- 
formation in an event to determine the most likely inter- 
action sequence. A photon of initial energy E (using the 
notation in Section 2) interacts in the instrument at N 
sites, depositing an energy of Ei at each location tv;. It is 



assumed that the interactions 1, . . . , N — 1 are Compton 
scatters, and interaction N is the final photoabsorption. 
Given the correct ordering of the interactions, there are 
two independent ways of measuring N — 2 of the scatter- 
ing angles, cos fa, ■ ■ ■ , cos 4>n-i ■ 

Geometrical measurement of cos fa. From simple vec- 
tor analysis, given the correct ordering of the interaction 
sites one can derive the scatter angles 



COS( 



= f i -f^_ 1 ,i = 2,...,N-l , 



(9) 



where the uncertainties in the scattering angles, S cos fa, 
can be estimated from the spatial uncertainty in the scat- 
tering angles (Equation 8), yielding 



5 (cos fa) = 5fa jT sin(/>j 



(10) 



Compton kinematics measurement of cos §\ . Given the 
correct ordering, the measured values of Wi can be de- 
rived, which were defined earlier as the energy of the pho- 
ton after each scattering i, in units of m e c 2 . The Compton 
scatter formula (Equation 2) gives: 
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Given the TV — 2 independent measurements of 
cos fa, . . . , cos0jv-i, a trial ordering of the interaction 
sites can be tested for consistency. If the assumed or- 
dering is incorrect cos fa will not equal cos fa in general. 
Every possible permutation of orderings can be tested to 
determine the one most consistent with cos^i = cosfai- 
Given a trial ordering, the two angle cosines for sites 
i = 2, . . . , N— 2 are relabelled for convenience rji = cos fa, 
Vi = cos fa,. 

As a first test, trial orderings that produce values of 
I rfi I— 1 are ruled out, since cos fa < 1 for any scatter- 
ing angle fa. This condition will eliminate many order- 
ings which cannot physically be due to multiple Compton 
scatters followed by photoabsorption. Next a least-squares 
statistic measuring the agreement between the redundant 
scatter angle measurements is defined: 



i N - 1 i 

i \ - (m - 

N - 2 S ^ 



Srtf) 



(13) 



In general, x 2 wm be minimized when the interac- 
tions are properly ordered (i.e. the order in which they oc- 
curred). Therefore, all possible permutations can be tested 
for their value of \ 2 , an d the ordering corresponding to the 
minimum value, Xmim ^ s taken as the most likely ordering. 

This consistency statistic also provides a powerful tool 
for rejecting background events. If the event is truly a 
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Fig. 4. Photopeak distributions for 3+ site events. The 
fraction of events with the first and second interaction 
sites spatially resolved are presented (o), along with the 
fraction of events which have been properly ordered (hence 
imaged) using CKD (A). The fraction of events rejected 
by the CKD statistic (□) as well as the fraction incorrectly 
imaged ( x ) are also shown. 



multiple Compton scatter event followed by a photoab- 
sorption then Xmin ~ 1- By setting a maximum accept- 
able level for Xmin> even ts that do not fit this scenario 
can be rejected. Such events include partially-deposited 
photons which scatter out of the instrument (Compton 
continuum) , photon interactions with spatially unresolved 
interaction sites, events with interactions below the de- 
tector threshold, pair-production events, and similarly /3 + 
decays. These events frequently have Xmin 1j allowing 
a strong rejection statistic that is not very sensitive on 
the level set on Xmin- Here, Xmin h as been treated as a 
normal least-squares statistic with N — 2 degrees of free- 
dom, and events are rejected which have probabilities of 
Xmin < 5%. Variations in the level between 1% and 10% 
do not strongly affect CKD rejection capabilities. For ex- 
ample, varying this level from 5% to 1% shifted the CKD 
efficiency curves in Figure 4 by 1-2%. 

The fraction of 3+ site photopeak events which have 
the first and second interaction sites spatially resolved - 
and hence could be imaged to the proper direction - is 
shown in Figure 4 as a function of energy, for the instru- 
ment model discussed in Appendix A. Roughly 90% of 
all events from 0.2-20 MeV have their first and second 
sites spatially resolved from each other. (Some of these 
events do not have their second, third, etc., interactions 
spatially resolved from each other, and will be rejected 
by the limits on Xmin-) This figure also shows the frac- 
tion of the photopeak events which CKD properly orders 
(correctly reconstructed), as well as the fractions improp- 
erly ordered (hence incorrectly imaged to off-source back- 
ground), and the fraction completely rejected. For ener- 
gies below ~10 MeV, CKD allows proper reconstruction 
(hence imaging) of ~ 60 — 70% of the photopeak events, 
while rejecting 10—20%. The remaining 10—20% are incor- 



rectly imaged into the off-source background. For compari- 
son, if the order of the interaction sites were randomly cho- 
sen < 15% would be correctly imaged, while the remaining 
> 85% would be incorrectly imaged into the background. 

5. Single Compton scatter events: single scatter 
discrimination 

CKD will only work for N > 2 since there are no indepen- 
dent scattering angle measurements for a single Compton 
scatter followed by a photoabsorption. It turns out that 
the ordering of two-site photopeak events can still be de- 
termined with a high probability; however, the ability to 
reject background events is lost. As is discussed further 
in Section 6, the loss of background rejection, coupled 
with low peak-to-Compton ratios, and a larger fraction 
of backscatter events mean that the inclusion of two-site 
events will likely hurt the sensitivity of a GCT; however, 
discussion of event ordering is still included for complete- 
ness. 

Given a two-site event, the first test one can perform 
is to determine whether both possible orderings of the in- 
teraction sites are energetically compatible with a single 
Compton scatter, i.e. are compatible with the requirement 
that cos(/»i < 1 (Equation 2). In Figure 5, the fraction 
of spatially-resolved, photopeak events with unique order- 
ings are plotted versus energy. Also plotted are the fraction 
with ambiguous orderings. At energies below ~0.4 MeV 
the majority of resolved photopeak events have a unique 
ordering, while at higher energies most events are ambigu- 
ous. 

As an empirical test of the ambiguous events, the rel- 
ative magnitude of the energy lost in the initial scatter 
(Ei) compared to the photoabsorption (E 2 ) can be com- 
pared. The fraction of resolved two-site photopeak events 
which have ambiguous orderings with Ei > E 2 is plotted 
in Figure 5. At higher energies, nearly all of the resolved 
photopeak events with ambiguous interaction orders have 
Ei > E 2 , which can be used to determine the most likely 
interaction order. This empirical result can be easily un- 
derstood, in hindsight, by the fact that photons which de- 
posit most of their energy in the initial Compton scatter 
are much more likely to be photoabsorbed in the second 
interaction. 

Therefore, a simple Single Scatter Discrimination 
(SSD) technique to determine the most likely interaction 
ordering of two-site events follows. First one determines 
whether a physically unique ordering exists; if not, the 
larger energy deposit is assumed to be the initial Compton 
scatter. Only at the lowest energies are two-site events pos- 
sible in which neither ordering is acceptable (unresolved 
events, Compton continuum), and some background re- 
jection is possible. 

The fraction of two-site photopeak events which are 
spatially resolved is shown in Figure 6 as a function of 
energy, for the instrument model in Appendix A. Roughly 



6 



Boggs et al.: Compton Event Reconstruction 



1 .0 - 

I 0.8 - 

o • 

o 

- o.6 : 

° 0.4 - 
c 
.9 

"u 0.2 1 

D 

0.0 j_ 
0.1 



□ A 



a A A g A a 



□ unique order 

a ambiguous order 

o E 1 >E 2 



□ □ 



1.0 10.0 
energy [MeV] 



6 r 

4 7 

2 7 

t_ 
0.1 



O spatially resolved - 
A correctly reconstructed 
□ rejected 

x improperly imaged 
<> 8 & * « " 



□ i □ i □ i i □ □□ □ □ ^ □ ^ □ □□ - 

1.0 10.0 
energy [MeV] 



Fig. 5. Empirical distributions for fully-resolved, two-site 
photopeak events. The fraction of events with only one 
physically possible ordering (□), as well as events with 
both orderings physically possible (A) are shown. Also 
shown are the events with both orderings possible, but 
with the larger energy deposit in the initial scatter rather 
than the photoabsorption, E\ > E$ (o). 



Fig. 6. Photopeak distributions for two-site events. The 
fraction of events with the first/second interaction sites 
spatially resolved are presented (o), along with the frac- 
tion of events which have been properly ordered (hence 
imaged) using SSD (A). The fraction of events rejected 
due to no physically possible orderings (□), as well as the 
fraction incorrectly imaged (x) are also shown. 



80% of all events from 0.2-20 MeV are resolved. This num- 
ber is about 10% lower than the 3+ site events, due to the 
smaller path lengths of the lower energy scattered pho- 
tons in two-site photopeak events. Also shown in Figure 6 
is the fraction of events which SSD has properly recon- 
structed. SSD allows proper reconstruction (hence imag- 
ing) of ~ 60 — 80% of the photopeak events, while improp- 
erly imaging the remaining ~ 20 — 40% into the off-source 
background. Only a relatively small number of low en- 
ergy events can be rejected outright. SSD is least effective 
around 0.5 MeV, where the unique/ambiguous ordering 
signatures are not as clear. For comparison, if the order of 
the interaction sites were randomly chosen ~ 40% would 
be properly imaged, while the remaining ~ 60% would be 
improperly imaged into the background. 



6. Full reconstruction with background rejection 

Given these two methods of determining the photon inter- 
action order in GCTs, CKD and SSD, a full event recon- 
struction technique incorporating other background rejec- 
tion techniques can be developed. The high spectral and 
spatial resolution of a GCT make several powerful back- 
ground rejection techniques possible. The predominance of 
multiple scattering events, while initially a complication, 
dramatically helps in the overall background rejection. 

The dominant sources of background in GCTs 
are expected to be diffuse cosmic 7-ray emis- 
sion, induced satellite 7-ray emission, and induced 
P + , f3~ radioactivities in the GeDs themselves 
Graham et al. 19971 pehrels 1985] , 



( Jean et al. 199 

Dean et al. 1991| , Naya et al. 1996 ). Source/background 
photons which scatter out of the instrument before 



depositing all of their energy, and hence are improperly 
imaged, must also be included in these calculations. 

Restrictions on the acceptable events can have a dra- 
matic effect on sensitivity of a GCT. Specifically, several 
factors can affect the angular resolution of the instrument 
as well as the background rates - such as the inclusion 
of backscatter events, limits on the accepted scatter an- 
gles, and the minimum acceptable lever arm - and must 
be included in any discussion of full reconstruction and 
background rejection. 

6.1. Shield veto 

Placing an active shield below the bottom GCT detector 
plane could be useful for rejecting background photons 
from below the instrument (induced satellite 7-ray emis- 
sion), as well as helping to reject Compton continuum and 
/3— decay events in the instrument. However, many of these 
events can be distinguished and rejected using CKD and 
other tests/restrictions outlined below; therefore, the use- 
fulness of including a shield in the GCT design must be 
studied in detail for a given telescope configuration. 

6.2. Restrictions on number of interaction sites: pair 
production/ [3 + decays 

While it is obvious that any single-site interactions should 
be rejected in a Compton telescope, events with ~8 or 
more interaction sites should also be rejected since these 
are very likely due to pair production events, as is evident 
from Figure 3, or similarly /3 + decays. (See Section 6.8) 
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6.3. CKD Xmin test: Compton continuum, unresolved 
interactions, etc. 



Using the tests outlined in Sections 4 and 5, the most likely 
ordering of the interaction sites can be determined, and for 
3+ site events many of the unresolved and Compton con- 
tinuum events, as well as pair production and (3— decays, 
rejected. 

Shown in Figure 7 is the peak-to-Compton ratio for 3+ 
site events, here defined as the ratio of the properly imaged 
photopeak events to the corresponding integrated Comp- 
ton continuum (photons which scatter out of the instru- 
ment before depositing all of their energy) . This standard 
measure for 7-ray spectroscopy instruments has an altered 
meaning here, since the Compton continuum events will 
be incorrectly imaged, and thus will appear as off-source 
background. The peak-to-Compton ratio is shown both 
before and after rejection of the continuum events with the 
CKD statistic. CKD rejection of the Compton continuum 
events increases the peak-to-Compton ratio by factors of 
3-6, an important improvement for low background instru- 
ments. By rejection of events appearing to originate from 
below the instrument (Section 6.4), as well as backscat- 
tered interactions (Section 6.5), this ratio can be increased 
by further factors of 2-4. 

Also shown in Figure 7 is the photopeak-to-Compton 
ratio for two-site events using SSD, which is significantly 
lower than the corresponding ratio for 3+ site events. 
In fact, this ratio drops significantly below unity, which 
means that more background than signal is being created 
in the instrument for two-site events. This result questions 
whether two-site events should be included in actual ob- 
servational analysis given the accompanying increase in 
unrejectable background. This conclusion is further sup- 
ported by the fact that the majority of two-site events are 
backscatters (Section 6.5), which will significantly degrade 
the angular resolution. Even though inclusion of two-site 
events is unlikely to improve the overall sensitivity for 
GCTs, detailed background analysis for specific instru- 
ment configurations is required to determine the overall 
effects. 



6.4. Effective TOF 

Once the most likely order of interactions and the initial 
scatter angle are determined, it is possible to determine 
whether the incident photon scattered upwards or down- 
wards in the instrument, as well as whether the initial 
scatter was forward or backwards. Thus, events which ap- 
pear to be photons originating from below the instrument 
can be rejected, which include the induced satellite 7-ray 
emission, many Compton continuum events which were 
not rejected by CKD, photons which scatter in the pas- 
sive satellite material before interacting in the detectors, 
and many of the pair production and (3— decay events. 
The simulation results for the configuration in Appendix A 
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Fig. 7. Peak-to-Compton ratios. Shown are the ratios of 
correctly ordered (imaged) photopeak events to integrated 
Compton continuum events for: 3+ sites before CKD re- 
jection (o), 3+ sites after CKD rejection (A), and two-site 
events (□). 

show that ~ 95% of photons originating from below the 
instrument are rejected. 

6.5. Backscatters 

Once the most likely ordering of interaction sites has 
been determined, this information can also be used to ac- 
cept/reject backscattered source photons. The fraction of 
photopeak events which backscatter during the initial in- 
teraction is not strongly energy dependent, ~ 60 — 70% 
for two site events, and ~ 30% for 3+ site events. These 
events can significantly increase the effective area at lower 
energies, where two-site events are most common, at the 
expense of degrading the angular resolution due to larger 
uncertainties in 8<pi,E for backscatters events (Equation 
4) . It is unlikely that the overall sensitivity will improve by 
including backscattered events given the increased back- 
ground rates and degraded angular resolution; however, 
the effects on sensitivity will depend on the exact instru- 
ment configuration and observational goals. 

6.6. 'Standard (f> restriction' 

Restrictions can also be set on the scattering angles ac- 
cepted for forward-scattering photons entering the front 
of the instrument. These limits can be used to re- 
strict the instrument FOV to improve imaging capabili- 
ties and background, such as the 'standard <f> restriction' 



(Schonfelder et al. 1982). These restricitions will have to 
be reanalyzed in detail for specific GCT configurations. 

6. 7. Nonlocalized (3~ decays 

In a nonlocalized (3~ decay, the daughter nuclide is 
produced in an excited state which quickly decays on 
timescales relative to the detector collection time, emit- 
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Table 1. Characteristic photon energies for the strongest 
nonlocalized (3~ decays in natural Ge. 



daughter 


photon energy 


75 Ge 


0.265 MeV 


73 Ga 


0.297 MeV 


72 Ga 


0.834 MeV 


71 Zn 


0.512 MeV 


76 As 


0.559 MeV 


28 Al 


1.779 MeV 


77 Ge 


0.216 MeV 



1 .0 -a 
0.8 7 
0.6 7 
0.4 7 
0.2 7 

0-0 L_ 
0.0" 



O0.5 MeV - 
A 2.0 MeV - 



7^ 

9h 



1.00 10.00 
ever arm [cm] 



100.00 



ting a photon with energy characteristic to the daughter 
nuclide. Therefore, the event consists of the intial (3~ de- 
cay site, plus the interaction sites of the emitted photon. 
Such an event can be rejected if the characteristic photon 
energy can be detected in any combination of the inter- 
action site energies. The seven dominant (3~ isotopes and 
characteristic photon energies for natural Ge are given in 
Table 1 jNaya et al. 1996| , |Gehrels 1985| ). In general, if the 
coincident 7-ray is fully deposited, then the event will have 
Xmin ~ 1> with the (3~ electron interaction ordered as the 
initial 'scatter' site. In these cases, W\ will have the char- 
acteristic photon energy specific to that decay. The rejec- 
tion of all events with W\ equal to one of the characteristic 
energies in Table 1 can dramatically decrease the (3~ decay 
background, with only a small effect, typically < 3% drop 
in photopeak efficiencies for true photon events. More (3~ 
decays can be rejected if every possible combination of 
interaction sites is tested for the decay photon energies, 
at the expense, however, of more rejection of photopeak 
events, typically ~ 15 — 20%. If the coincident 7-ray is 
only partially deposited, the event will likely be rejected 
by the limits set on Xmin- 

f}~ decay background events were simulated for the 
instrument discussed in Appendix A. The results of these 
background calculations will be presented in a separate pa- 
per, but here we make preliminary use of these simulations 
to demonstrate the background rejection capabilities. Af- 
ter initial rejection of two site (34.4%) and 8+ site (2.9%) 
events, 62.7% of the nonlocalized (3~ events remain. Ap- 
plying the CKD test, requiring a 5% probability of Xmim 
brings the remaining number of decays down to 17.9%. Af- 
ter screening the interactions for characteristic /3~ decay 
energies, this number is reduced to 9.3%(15.1%). Finally, 
after rejecting events which appear to originate from be- 
low the instrument, or which appear to be backscatter 
events, the final number of unrejected (3~ decays comes 
to 4.2%(6.8%), a factor of 20 (15) reduction in this back- 
ground component. (First numbers give results when all 
combinations of interaction sites are searched for (3~ de- 
cay energies, while the numbers in parenthesis are results 
when only W\ is tested. Typical errors ~ 0.2%.) 



Fig. 8. The fraction of 0.5 and 2.0 MeV photopeak events 
with lever arms (separations between the first and second 
interaction sites) greater than the specified values for the 
model in Appendix A. Shown for comparison are several 
characteristic distances of this model. 



6.8. Positron signatures 

A further test for rejecting the pair production//3 + back- 
ground events that survive the other tests/restrictions 
outlined above is to search for positron annihilation sig- 
natures in the interaction energies. By analyzing all com- 
binations of the interaction sites to see if the energies sum 
to m e c 2 = 0.511 MeV, events with a positron annihilation 
signature can be rejected. This test typically reduces the 
non-pair production photopeak events by < 2%. 

P + background events were simulated for the instru- 
ment discussed in Appendix A. After initial rejection of 
two site events (23.0%) and 8+ site events (3.0%), 74.0% 
of the events remain. After the CKD test, 16.7% remain. 
After screening the interactions for 0.511 MeV positron 
annihilation signatures, the number is reduced to 5.2%. Fi- 
nally, after rejecting events which appear to originate from 
below the instrument, or which appear to be backscatter 
events, the final number of unrejected (3 + events is 1.9%, 
a factor of 50 reduction in this background component. 
Similar reductions occur when these tests are applied to 
pair production events. (Typical errors ~ 0.1%.) 

6.9. Minimum lever arm 

In general, a minimum acceptable distance between the 
first and second interaction sites - the lever arm - must 
be set. Figure 8 shows the fraction of 0.5 and 2.0 MeV 
photopeak events with lever arms above a given level, for 
the instrument configuration in Appendix A. Similar to 
the case of backscattered events, a smaller minimum lever 
arm means a higher effective area at the expense of poorer 
angular resolution. The exact lever arm chosen will depend 
on the instrument configuration and observational goals. 
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Table 2. Percentage of events remaining after subsequent application of rejection techniques. 
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rejection 


0.5MeV 


2.0MeV 


8 decays (3 


+ decays spacecraft 


technique 


photopeak photopeak 








2 site events 


65.0% 


82.0% 


65.6% 


77.0% 


61.6% 


8+ site events 


64.9% 


80.9% 


62.7% 


74.0% 


60.8% 


CKD 


51.2% 


61.4% 


17.9% 


16.7% 


35.2% 


P signatures 


48.4% 


58.0% 


15.1% 


5.2% 


32.2% 


backscatter/TOF 


35.5% 


38.5% 


6.8% 


1.9% 


12.2% 


min lever arm (10 cm) 


15.0% 


16.7% 


2.0% 


1.2% 


2.6% 



7. Conclusions 

Event reconstruction in future high resolution Compton 
telescopes will present a number of complications com- 
pared to historical configurations. The initial complication 
of multiple-scattering of photons in GCTs, however, turns 
out to be an advantage: the application of CKD to 3+ 
site events, combined with the high spectral and spatial 
resolution of GcDs, allows extremely efficient background 
suppression, crucial for Compton telescope performance. 
This paper has outlined a set of tests and restrictions, ac- 
counting for realistic instrument/detector performance, to 
reconstruct photopeak events in GCTs while rejecting a 
large fraction of the background events. Table 2 presents 
the fraction of events, photon and background, that re- 
main after each rejection technique is subsequently ap- 
plied. (The numbers in Table 2 assume only W\ is tested 
for (3~ decay energies.) Development of these event recon- 
struction techniques allows realistic evaluation of the per- 
formace and sensitivity of GCT designs. Our next goal is 
to simulate the efficiency, resolution, background and sen- 
sitivity of several Compton telescope configurations, uti- 
lizing the event reconstruction techniques developed here 
to realistically determine the performace of these instru- 
ments. CKD rejection has been shown to be the most ef- 
ficient background rejection technique; however, the addi- 
tion of effective TOF, backscatter, nonlocalized (3~ decay, 
and positron signature tests dramatically improve back- 
ground rejection capabilities. We anticipate that use of 
these techniques will achieve overall sensitivity improve- 
ments in GCTs by factors of ~ 5 — 10. 
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Appendix A: Example GCT configuration 

While it is not the intention of this paper to fully char- 
acterize the performance of a specific GCT, it is useful 
to have a telescope model for which the results of the 
event reconstruction can be presented. The telescope con- 
figuration modeled in this study is presented in Figure 1. 
The instrument consists of five planar arrays of 15 mm 
thick germanium, each of area 100 cm x 100 cm. In real- 
ity each array would consist of separate smaller detectors 
(<~ 5 cm x 5 cm) tiled to form the entire plane; however, 
the simulation performed here modeled each plane as a 
solid detector for simplicity. The five planar arrays are 
spaced 20 cm apart. 

This configuration differs from historical Compton 
telescope configurations which generally consist of two de- 
tector planes separated by 100-150 cm. This separation 
distance is determined by the spatial resolution in z and 
the desired angular resolution. As will be discussed in a 
second paper, the configuration modeled here significantly 
improves the effective area of the telescope by letting each 
plane act as converter, and permitting a much wider range 
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Fig. A.l. Assumed resolutions for the GCT model. (Left) 
The total energy resolution (solid line) is determined by 
the electronic noise (dotted line) and the intrinsic resolu- 
tion (dashed line) added in quadrature. (Right) The total 
x-, y- (solid line) and z- (dot-dashed line) spatial resolu- 
tion are determined by the segmentation/signal processing 
limits (dotted line), and the recoil electron range (dashed 
line) added in quadrature. (See text for details.) 



of scatter angles to produce good events. Allowing large- 
angle scatters also significantly increases the instrument 
FOV, and limits the effects of point spread function smear- 
ing for sources at large off-axis angles. The potential draw- 
backs of this configuration are increased background and 
degraded angular resolution. 

The instrument was simulated using CERN's GEANT 
Monte Carlo code. The Monte Carlo simulation produces 
a file of interaction locations and energy depositions for 
each photon//?— decay event. Before performing event re- 
construction on the interactions, the simulated events are 
modified to reflect realistic measurement uncertainties of 
an instrument: for each interaction, a random Gaussian- 
distributed uncertainty is added to the energy and posi- 
tion of each interaction. All interaction locations which lie 
within twice the instrumental spatial resolution of each 
other are combined into a single interaction site, to accu- 
rately reflect the resolving power of the detectors. Finally, 
interaction sites with energy deposits below the assumed 
detector threshold of 10 keV are ignored. 

Two components are assumed to add in quadrature 
to determine the energy resolution: (i) a constant elec- 
tronic noise, W e — 1.0 keV FWHM, and (ii) the intrin- 
sic resolution Wi determined by the germanium Fano fac- 
tor, F = 0.13, and average free electron-hole pair energy, 
e = 2.98 eV, giving Wi = 2.35V¥eE FWHM. This cor- 
responds to a resolution ~ 1.8 keV FWHM at 1 MeV, 
which is optimistic but not unrealistic. It is assumed that 
charge trapping and ballistic deficit do not significantly 
alter this energy resolution. The two components as well 
as the total energy resolution are shown in Figure Al. 

It is assumed that two components add in quadra- 
ture to determine the 1 — D spatial resolutions, Sx, Sy, 5z, 



of the detectors: (i) the range of the recoil electrons in 
the detector, and (ii) the positioning limits of the de- 
tector due to physical segmentation and/or signal anal- 
ysis. Calculated electron ranges in germanium for differ- 



ent energies (Mukoyama 1976) are used here as the 1- 
D FWHM positional uncertainties, Sx e ,Sy e ,5z e . Meth- 
ods to determine the event position by physically seg- 
menting the GeD contacts into cross strips or pixels 



(Luke et al. 1994, Kroeger et al. 1996), as well as using 



advanced signal processing to interpolate to even better 
positions ( Boggs 199g , Luke et al. 1994 ), are currently ac- 
tive fields of research - so this component of the spatial 
resolution remains speculative for now. Here it is assumed 
that signal processing will allow positional resolutions of 
~ 0.5 mm FWHM at 100 keV, and that the discrimi- 
nation capabilities go as the signal-to-noise ratio of the 
induced detector signal to electronic noise, i.e. as the in- 
verse power of the interaction energy. It is also assumed 
that there is ~ 1 mm physical segmentation of the de- 
tector contacts in x, y, so that this component never ex- 
ceeds this value. The z uncertainty, however, is not con- 
strained by any such segmentation at the lowest energies. 
Therefore, the signal processing uncertainty is given by 
Sx Sl 5y s ,Sz s ~ 0.50(£:/100/fceV r )- 1 mm FWHM, maxi- 
mizing at 1 mm in x, y below 50 keV, and approaching, 
but never maximizing at 15 mm in z at low energies. The 
two components as well as the total spatial resolution are 
shown in Figure Al. 



